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ABSTRACT 

PSRJ1811-1736 (P s = 104ms) is an old (~ 1.89 Gyrs) binary pulsar (P orb =18.8 d) in a 
highly eccentric orbit (e = 0.828) with an unidentified companion. Interestingly enough, the 
pulsar timing solution yields an estimated companion mass 0.93 M® ^ Mc ^ 1.5 Mq, 
compatible with that of a neutron star. As such, it is possible that PSR J181 1 — 1736 is a 
double neutron star (DNS) system, one of the very few discovered so far. This scenario can be 
investigated through deep optical/infrared (IR) observations. We used J, H, A'-band images, 
obtained as part of the UK Infrared Telescope (UKIRT) Infrared Deep Sky Survey (UKIDSS), 
and available in the recent Data Release 9 Plus, to search for its undetected companion of the 
PSR J181 1 — 1736 binary pulsar. We detected a possible companion star to PSR J181 1 — 1736 
within the 3cr radio position uncertainty (1'.'32), with magnitudes J=18.61 ± 0.07, H=16.65 ± 
0.03, and K= 15.46 ± 0.02. The star colours are consistent with either a main sequence (MS) 
star close to the turn-off or a lower red giant branch (RGB) star, at a pulsar distance of ~ 5.5 
kpc and with a reddening of E(B — V) « 4.9. The star mass and radius would be compatible 
with the constraints on the masses and orbital inclination of the binary system inferred from 
the mass function and the lack of radio eclipses near superior conjunction. Thus, it is possible 
that it is the companion to PSR J181 1 — 1736. However, based on the star density in the field, 
we estimated a quite large chance coincidence probability of ~ 0.27 between the pulsar and 
the star, which makes the association unlikely. No other star is detected within the 3er pulsar 
radio position down to 20.5, H^ 19.4 and 18.6, which would allow us to rule out a 
MS companion star earlier than a mid-to-late M spectral type. 
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1 INTRODUCTION 

The radio pulsar PSRJ1811-1736 (P a = 104 ms) was detected at 
1374 MHz (Lyne et al. 2000) during the Parkes Multibeam Pulsar 
survey (Manchester et al. 2001). It is in a binary system, with an 
orbital period P or b=18.8d and a high eccentricity e=0.828 (Coro- 
ngiu et al. 2007). The updated timing parameters, including gen- 
eral relativistic effects, give a period derivative P s ~ 9.01(5) x 
10~ 19 ss _1 which yields a spin-down age tsd ~ 1-89 x 10 9 yr 
and a surface magnetic field B sur f ~ 9.8 x 10 9 G. The P B and P s 
suggest that PSR J181 1 — 1736 is a mildly-recycled pulsar, i.e. the 
spin-up phase via matter accretion from the companion star was 
too short for the pulsar to reach a spin period of a few ms, typi- 
cal of fully-recycled pulsars. A possible scenario is that the com- 
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panion was an high-mass star which underwent a supernova ex- 
plosion, itself turning into a neutron star (Bhattacharya & van den 
Heuvel 1991). Thus, PSRJ1811-1736 might be one of the ss 10 
double neutron star (DNS) systems out of the w 2000 radio pulsars 
known to date (Manchester et al. 2005). The DNS picture is rein- 
forced by the limits on the companion mass, derived from the mass 
of the system Mtot = 2.57 ± 0.10 A/© inferred from the mea- 
surement of the periastron advance lj — 0.0090° ± 0.0002° yr -1 
and from the mass function. For a pulsar mass Mp ^ 1.17 M®, 
larger than the minimum value inferred for a radio pulsar (PSR 
J15 18+4904; Janssen et al. 2008), this yields a companion mass of 
0.93 M© ^ Mc < 1.5M© (Corongiu et al. 2007), compatible 
with that of a neutron star. 

A conclusive piece of evidence that PSR J181 1 — 1736 is a 
DNS would be the detection of its companion as a radio pulsar, 
like in the double pulsar PSR J0737-3039A/B (Lyne et al. 2004). 
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However, it escaped detection so far, perhaps because of an un- 
favourable beaming or because it is no longer in its active radio 
phase. Alternatively, a conclusive piece of evidence would be the 
non-detection of the companion in deep optical/infrared (IR) ob- 
servations. The pulsar companion is not detected in the Digitised 
Sky Survey (DSS) down to R » 22 (Mignani 2000) and in the 2 
Micron All Sky Survey (2MASS) down to K s « 15, computed at 
the Lyne et al. (2000) and Corongiu et al. (2007) radio positions, 
respectively, with the latter limit being quite uncertain owing to the 
much higher crowding in the pulsar field at IR wavelengths. Such 
limits would only rule out a giant companion but, for the allowed 
mass range they would still be compatible with a mid to late-type 
main sequence (MS) star, a white dwarf, or a neutron star. No deep 
optical/near-IR observations of PSR J 18 1 1 — 1736 have ever been 
performed so far. As suggested in Mignani (2000), given the sub- 
stantial interstellar extinction towards the pulsar near-IR observa- 
tions are more suited than the optical ones to set constraints on the 
companion star. 

Here, we present the results of a new investigation of the PSR 
J181 1 — 1736 field using IR survey data much deeper than 2MASS. 
The observations and results are discussed in Sectn. 2, while the 
implications for the PSR J181 1 — 1736 companion are discussed in 
Sectn. 3. 




2 INFRARED OBSERVATIONS AND RESULTS 
2.1 Observation description 

No near-IR observations of the PSRJ1811 — 1736 are available in 
either the EScQ or the Geminjf] Science Data Archives. Thus, we 
searched for near-IR data of the PSR J 1 8 1 1 — 1736 field in the im- 
age archive of the UK Infrared Deep Sky Survey (UKIDSS) per- 
formed with the Wide Field Camera (WFCAM) at the UK Infrared 
Telescope (UKIRT) at the Mauna Kea Observatory (Hawaii). WF- 
CAM (Casali et al. 2007) is a mosaic detector of four 2048x2048 
pixel Rockwell devices, with a pixel scale of 0'.'4 and covering a 
field-of-view of 0.21 square degrees. A general description of the 
UKIDSS survey is given in Lawrence et al. (2007). The UKIDSS 
survey covers several regions, with a different sky coverage, and 
sensitivity limits in the ZYJHK UKIRT photometric system 
(Hewett et al. 2006). The field of PSR J1811 1736 is included 
in the Galactic Plane Survey (GPS; Lucas et al. 2008) which cov- 
ers about 1800 square degrees in JHK down to sensitivity limits 
which are more than a factor of ten deeper than 2MASS. Like all the 
UKIDSS data, the GPS images are processed through a dedicated 
pipeline (Hambly et al. 2008) developed and operated at Cam- 
bridge Astronomical Survey Unit (CASU) which performs basic 
reduction steps (dark subtraction, flat fielding), image de-jittering, 
stacking, and mosaicing. The pipeline also runs a source detection 
algorithm and produces source catalogues. Astrometry and pho- 
tometry calibration are performed using 2MASS stars as a refer- 
ence (Hodgkin et al. 2009). We searched for the reduced science 
images of the PSR J 18 1 1 — 1736 field and associated object cat- 
alogues through the WFCAM Science Archive (WSAfj interface 
accessible via the Royal Observatory Edinburgh. We queried the 
most recent UKIDSS Data Release (version 9 Plus) made available 
on October 25th 2011. The field was observed on July 18 2006. 
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Figure 1. (top) V X 1' K s -band image of the PSR J1811-1736 field ob- 
tained from 2MASS. The 20" X 20" square corresponds to the sky area 
shown in the bottom panel, (bottom) K s -band image zoom of the same field 
obtained from UKIDSS. North to the top, East to the left. The radio po- 
sitions of PSR J181 1 — 1736 computed with and without fitting its proper 
motion in the timing model (see Sectn. 2.2) are marked by the two sets of 
circles, drawn with thick and thin lines, respectively. In both cases, the inner 
circles correspond to the lcr (()'.' 44; 0" 19) uncertainty radii, while the outer 
circles correspond to the 3cr (l'.'32; 0''58) ones. The star within the proper 
motion-corrected 3cr radio error circle (unresolved in the 2MASS image) 
has magnitudes J=18.61±0.07, H=16.65±0.03, and K s = 15.46±0.02. 



We downloaded 10' x 10' J, H, if -band stacks around the pulsar 
position and the associated multi-band object catalogues. 



2.2 Pulsar astrometry 

For the search for the companion star to PSR J181 1 — 1736, we as- 
sumed as a reference its radio timing coordinates. We note that the 
pulsar radio timing solution presented by Corongiu et al. (2007) 
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Figure 2. K-(H-K) (left) colour-magnitude diagrams of the PSR J181 1-1736 field obtained from the UKIDSS photometry of stars detected in a 10' X 10' 
region around the PSR J181 1 — 1736 position. The star detected closest to the pulsar position (see Fig. 1) is marked in red. Average photometric errors in the 
pulsar field, for different magnitude and colour bins, have are also plotted in the left and right panels respectively. (K — H) vs. (J — K) colour-colour diagram 
of the same region. The reddening vectors for an E(B — V) = 3 are shown, as a reference. 




Figure 3. Same as Fig. 2 but with UKIDSS data for the Baade's Window (light grey) overlaid for comparison. 



is based on data taken in the epoch range MJD=50842-53624 and 
does not include the determination of the pulsar's proper motion, an 
essential parameter for recomputing its position at a given epoch. 
For this reason, we re-analysed the data presented in Corongiu et 
al. (2007) adding the proper motion to the timing model, to obtain 
a new radio timing position at a reference epoch (MJD=53624) 
closest to that of the UKIDSS observation (MJD= 53934). Thus, 
we obtained aj20Q0 = 18 h ll m 55=0385 ± 0:0029 and feooo = 
-17°36'38'.'45±0'.'41forfhe position and ^ a cos(S) = 18.3±9.7 
masyr -1 and fj,$ = —176 ± 100 masyr -1 for the proper motion. 



where all quoted uncertainties are at lcr level. The extrapolated 
timing position at the epoch of the UKIDSS observation (MID= 
53934) is, then: a?j2ooo = 18 h ll m 55=054 ± 0!009 and S 

J2000 = 

-17° 36' 38'.' 60 ± 0'.'42, with an uncertainty radius of 0'.'44 (lcr) 
that accounts for the position uncertainty due to the proper motion 
extrapolation. For comparison, by applying the same timing model 
as above but without adding the proper motion, we obtain aj2ooo = 
18 h ll m 55=0337 ± 0:0014 and <5 J2 ooo = -17° 36' 37'.' 80 ± 0'.'19, 
where the choice of the reference epoch (MJD=53624) within the 
range spanned by the timing data is, in this case, arbitrary. Although 
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this position has a nominal uncertainty radius (0'.'19; la) that is 
smaller than obtained in the previous case, assuming it as a refer- 
ence at the epoch of the UKIDSS observations would introduce an 
unknown systematic uncertainty due to the neglected pulsar proper 
motion. For this reason, it is formally less correct than the radio 
position obtained by fitting the proper motion, despite the latter 
having a larger uncertainty radius. Nonetheless, in the following 
section we conservatively consider both positions in our search for 
the PSR J181 1 — 1736 counterpart. In computing the overall uncer- 
tainty on the PSR J1811 1736 position in the UKIDSS images 
we also accounted for systematics associated with the nominal ac- 
curacy on the UKIDSS astrometry calibration (0'.'05 rms at low 
Galactic latitudes; Lawrence et al. 2007), the internal astrometric 
accuracy of 2MASS (<072 for stars with 15.5 < K «S 13), and 
the accuracy on the link of 2MASS to the International Celestial 
Reference System (f/.'OlS; Skrutskie et al. 2006). 

2.3 Results 

The UKIDSS A'-band image of the PSR 11811-1736 field is 
shown in Fig. 1 (bottom) compared to the corresponding 2MASS 
image (top). For comparison, we plotted the two pulsar positions 
derived from the radio timing solution, with and without fitting the 
proper motion. As seen, no object is detected within the two la ra- 
dio position error circles (0'.'44 and 0'.'19 radii, respectively). How- 
ever, a star (K = 15.46 ± 0.02), unresolved in the 2MASS image 
but clearly detected in the much higher resolution UKIDSS one, 
is detected within the proper motion-corrected 3a radio error cir- 
cle (1'.'32 radius). Thus, its association with the pulsar cannot be 
ruled out a priori and needs to be investigated. The star is also de- 
tected in the J and H bands, with magnitudes of J = 18.61 ± 0.37 
and H — 16.65 ± 0.03. No other star is detected at, or close to, 
the computed 3a radio pulsar positions down to 3a limiting mag- 
nitudes of J~ 20.5, H~ 19.4 and K~ 18.6, as computed from 
the rms of the sky background (Newberry 1991). Given the high 
star density along the Galactic plane, however, the match can be 
the result of a chance coincidence. We computed this probability 
as P = 1 — exp(— ivpr 2 ), where r (~ l'.'32) is the matching ra- 
dius, assumed equal to the 3a uncertainty on the proper motion- 
corrected pulsar radio position, and p is the density of stellar ob- 
jects within an area of 10' x 10' around the pulsar. We found that 
p ~ 0.057 arcsec -2 , which gives P ~ 0.27. Such a high chance 
coincidence probability suggests that the star is likely unrelated to 
the pulsar, although we need a direct piece of evidence to firmly 
rule out the association. 



3 DISCUSSION 

3.1 The interstellar extinction in the pulsar direction 

We investigated whether the characteristics of the star detected 
close to the PSR J181 1 — 1736 position are compatible with it be- 
ing its companion star. To this end, we tried to determine its spectral 
type from its colours. Fig. 2 shows the colour-magnitude diagram 
(CMD) (H — K)-K and the colour-colour diagram (K - H) vs. 
(J — K) built from the photometry of field stars detected within 
a 10' x 10' area around the pulsar position, as derived from the 
UKIDSS object catalogue. We also plotted the location of the star 
detected close to the PSR J1811-1736 position (Fig. 1, bottom), 
whose location in both diagrams is consistent with the sequence of 
field stars. Thus, determining its spectral type from the comparison 



of its colours and flux with those of field stars is not obvious. More- 
over, the determination of the star's intrinsic colours is affected by 
the substantial interstellar extinction towards the pulsar, which is 
located in the Galactic plane (/ = 12.828°; 6 = 0.435°). In par- 
ticular, the CMD is very broadened, suggesting that the field is af- 
fected by a quite high, and probably differential, extinction. 

The interstellar extinction towards PSR J 1 8 1 1 — 1736 is uncer- 
tain, and this affects our estimate of the upper limits on the com- 
panion star luminosity. A first estimate of the interstellar extinction 
can be derived from the integrated Hydrogen column density along 
the line of sight to the pulsar. This is 7Y H = (1.24-1.64) x 10 22 , as 
computed using the HEASARC tool webpimm^l according to the 
Dickey & Lockman (1990) and Kalberla et al. (2005) Hydrogen 
maps. This gives E(B — V) = 2.2-2.9, according to the relation 
of Predhel & Schmitt (1995). However, PSR J181 1-1736 is closer 
than the edge of the Galaxy, at a distance D — 5.71q'^ kpc, es- 
timated from the radio pulse dispersion measure (DM=476 ± 5 pc 
cm" ; Corongiu et al. 2007) and the Galactic free electron density 
along the line of sight (Cordes & Lazio 2002). This would suggest 
a lower interstellar extinction. According to the Galactic extinction 
maps of Hakkila et al. (1997), the pulsar distance and Galactic co- 
ordinates would imply an interstellar extinction E(B — V) ~ 1.9. 
However, these estimates are only indicative mainly because of 
the uncertainties on the extinction maps on smaller angular scales. 
Unfortunately, the PSR J181 1 — 1736 field has not been observed 
in X-rays, so that no independent measurement of the interstel- 
lar extinction can be inferred from the hydrogen column density 
Nn directly derived from the fits to the X-ray spectra. In princi- 
ple, an independent measurement of the JVh can be obtained from 
the DM itself assuming an average ionisation fraction of the in- 
terstellar medium (ISM) along the line of sight. In the case of 
PSR J1811-1736, a DM=476 ± 5 pc cm" 3 would correspond to 
Nu ~ 1.5 x 10 22 cm -2 , for a 10% ionisation fraction. This would 
imply an E(B — V) ~ 2.7. However, this method is usually applied 
to pulsars closer than ~ 300 pc (e.g., Pavlov et al. 2009; Tiengo et 
al. 201 1) and is intrinsically affected by a much larger uncertainty 
for pulsars at larger distances, such as PSR Jl 81 1 — 1736. 

We tried to derive an independent estimate on the redden- 
ing along the line of sight by comparing the CMDs and colour- 
colour diagrams of field stars with those in a reference region of 
very low reddening, such as the Baade's Window. As we did for 
the PSR J1811-1736 field, we extracted from the UKIDSS data 
the object catalogues relative to a 10' x 10' area centred around 
the Baade's Window, for which we assumed the coordinates of 
the globular cluster NGC 6522: a. J2 ooo = 18 h 03 m 34 ! !08 and 
feooo = -30°02'02'.'3 (Di Criscienzo et al. 2006). Fig. 3, shows 
the same diagrams as in Fig. 2 but with the UKIDSS data for the 
Baades Window region overlaid. From the comparison of the two 
sets of diagrams, we derived an estimate of the interstellar extinc- 
tion towards the pulsar. Firstly, we computed the average of the 
distribution in the colour-colour space for the Baade's Window re- 
gion, applying a 3a clipping. Secondly, we did the same for the 
pulsar field but selecting a region of 50" radius around the pulsar 
position not to be affected by the differential extinction in the field. 
Then, from the difference between the two values we estimated an 
E(B — V) — 5.7 ± 1.9 for the pulsar field. This value is a factor 
of 2 larger than the Galactic interstellar extinction inferred from the 
Hydrogen column density maps in the pulsar's direction. However, 
we note that the Nn value reported above is a weighted average rel- 
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Figure 4. Measure of the granularity on the PSR J181 1 — 1736 field (y pa- 
rameter) computed using the method described in Gosling et al. (2006) as 
a function of the angular distance from the pulsar. The three curves cor- 
respond to the granularity measured in the J, H, and K s -band images (see 
legenda). 

ative to a 1° radius area around the pulsar coordinates, which does 
not rule out the presence of patches of higher Hydrogen column 
density on angular scales smaller than 0.4°, which are not resolved 
by the available maps. 

Indeed, it has been found that the interstellar extinction to- 
wards the Galactic bulge region is not uniform and shows strong 
variations, or granularities, on angular scales as small as 1' (see, e.g. 
Gosling et al. 2006). We used the UKIDSS J, H, and K s -band im- 
ages of the PSR J181 1 — 1736 field to measure the granularity of the 
interstellar extinction in the region, following the method described 
in Gosling et al. (2006). We considered a region of 500" x 500" 
centred on the PSR J181 1 — 1736 position. We divided the region 
in cells with dimension variable between 8" and 500" to sample 
the granularity in the field on different angular scales. For each cell 
dimension, we calculated the y parameter (see Eqn. 1 in Gosling et 
al. 2006), which gives a quantitative estimate of the granularity of 
the field and is defined as the variance of the number of stars in all 
cell normalised to the mean number of stars per cell. We computed 
the y parameter for the J, H, and K a -band images. For the region 
considered in our analysis, we found that the y parameter, hence 
the granularity, decreases as a function of wavelength (see Fig. 4), 
as in the case of one of the test fields used by Gosling et al. (2006). 
Thus, the high level of granularity in the pulsar field seems to be 
correlated with a high and variable reddening, explaining the large 
scatter in the CMD and colour-colour diagram of the field stars 
(Fig. 2). We note that the measured angular scale of the granular- 
ity in the PSR Jl 8 1 1 — 1736 field (see Fig. 4) is comparable to the 
50" radius region that we used to estimate the extinction in the pul- 
sar's direction from the CMD and colour-colour diagram analysis 
(see previous paragraph). Thus, we are confident that our procedure 
does not under/overestimate the assumed extinction value along the 
line of sight to the pulsar. 

3.2 The candidate companion star 

Under the hypothesis that the star seen close to the PSR 
J181 1 — 1736 position is its companion, we tried to determine its 
spectral type assuming the range of reddening values computed 
above. We considered a range of distances within the computed la 
uncertainty range on the PSR J181 1 — 1736 distance based on the 
DM (15=4.9-6.54 kpc), and an age range of 1-13.2 Gyr, consistent 
with the ages of the stellar populations in the Galactic centre re- 
gion (e.g. Zoccali et al. 2003), where the pulsar is located. We also 



considered different values of the metallicity Z. For older popula- 
tions, we considered both Z = 0.008 and Z — 0.02, while for 
the younger populations we considered only Z = 0.02. Then, for 
different values of metallicity, age, and distance we determined the 
best E(B - V) values in the range E(B - V) = 3.8-7.6 that 
minimise the sum of the projected distances of the star location 
in the J-( J — H) and K-(H — K) CMDs from the isochrones 
computed from the stellar models of (Marigo et al. 2008; Girardi et 
al. 2010). From the best combinations of metallicity, age, distance, 
and reddening, we then derived the corresponding mass (Mc) and 
radius (Rc) of the candidate companion star to PSR J181 1 — 1736, 
from the comparison with the isochrones. Finally, we combined 
the computed mass of the candidate companion star with the mass 
function of the pulsar system (Corongiu et al. 2007) and its total 
mass (M to t = 2.57 ± 0.10 Mq). In this way, we derived the pul- 
sar mass (Mp) and the orbital inclination angle i of the system for 
each combination of metallicity, age, distance, reddening, and mass 
of the candidate companion star. We filtered out combinations for 
which it resulted sini > 1. We also selected the acceptable com- 
binations for the conditions that the companion mass 0.93 Mq SC 
Mc < 1.5M and the pulsar mass 1.17 M Q < M P < 1.6 Mq, 
as it results from the PSR J 1 8 1 1 — 1736 timing analysis (Corongiu 
et al. 2007), where the lower limit on Mp corresponds to the min- 
imum measured value for the mass of a neutron star (Janssen et al. 
2008). The different parameter combinations are summarised in the 
first eight columns of Table 1 and 2. These parameters are consis- 
tent, for a reddening E(B — V) = 4.5-5, with a companion star 
still on the main sequence (MS) and close to the turn-off or on the 
lower Red Giant Branch (RGB), with an inferred mass Mc ~ 1— 
1.3 M e and radius R c « 4.3-5.8 R e . 

All the combination of parameters reported in Tab. [T] and |2] 
have been further checked against the lack of eclipses at superior 
conjunction (orbital phase (f> = 0.25) in the radio timing observa- 
tions. This check is grounded on the fact that the pulsar cannot be 
eclipsed at a given orbital phase if the corresponding pulse's time of 
arrival (ToA) has been determined, since ToA determination strictly 
requires the detection of the pulse. Hence, we calculated the orbital 
phases for each ToA presented in Corongiu et al. (2007), and we 
obtained that the closest available ToAs before and after superior 
conjunction correspond to an orbital phase (f> = 0.191366 and <j> 
= 0.277151, respectively. A parameters' combination is acceptable 
only if, at both orbital phases, the projected distance d of the pulsar 
to the centre of the companion, computed on the plane perpendic- 
ular to the line of sight, is larger than the companion radius, i.e. 
d/Rc > 1. Columns 9 and 10 of Tab.Q] and [2] report the values 
of d/Rc for the two values of the orbital phase computed above, 
with the possible combinations flagged YES and NO for the cases 
d/Rc > 1 and d/Rc < 1, respectively. Our calculation shows 
that for only 4 out of the 63 possible combinations of selected pa- 
rameters the required condition is satisfied at both orbital phases. 
These combinations imply a ~ 5 Gyr companion star, with mass 
Mc « 1.3Mq and radius Rc ~ 5-Rq, at a distance of ~ 5.5 
kpc and with a reddening E(B — V) » 4.9. This corresponds to a 
pulsar mass Mp w 1.3Af and inclination angle for the system of 
« 45° . Thus, according to the constraints on the masses and orbital 
inclination of the binary system, it is theoretically possible that the 
star detected at the radio position is, indeed, the companion to the 
pulsar. In this case, PSR J1811 — 1736 would not be aDNS. 

A ~ 5 Gyr MS companion would be compatible with the 
pulsar spin-down age (tsd = 1.9 Gyr) but not much so with the 
recycling scenario and the pulsar orbital parameters. In principle, 
the high orbital eccentricity of PSR J181 1-1736 (e = 0.828) can 
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Table 1. Best parameter combinations for the candidate companion star of PSR Jl 81 1 — 1736 obtained from the comparison between its location in the CMDs 
and model isochrones. Columns list the star metallicity Z, age (in logarithmic units), distance (D), reddening, its mass (Mq) and radius (Rq), the mass of the 
pulsar (Mp) in solar units, and the system inclination angle (i), obtained from the estimated companion mass, the system mass function, and the total mass. 
The next two columns list the projected distance (d) of the pulsar to the centre of the companion on the plane perpendicular to the line of sight, normalised to 
the companion radius Rq for the two values of the orbital phase <j> when the pulsar is observed closest to the superior conjunction <f> = 0.25. The last column 
flags the acceptable configurations under the condition that d/Rc > 1 at both orbital phases (see text). 
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4.99 
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NO 


0.020 
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NO 


0.020 
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4.41 
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NO 


0.020 
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4.71 
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1.54 
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0.81169 


0.55665 


NO 


0.020 
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4.75 


1.03 


5.03 


1.54 
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0.92585 


0.63339 


NO 


0.020 
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4.99 


4.78 


1.03 


4.41 


1.54 
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1.05601 


0.72244 


NO 
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5.78 
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0.58041 


NO 


0.020 
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5.70 


4.75 


1.05 


5.03 


1.52 


64.72 


0.95635 


0.66695 


NO 


0.020 


10.00 


4.99 


4.79 


1.05 


4.42 


1.52 


64.72 


1.08834 


0.75900 


NO 


0.008 


10.00 


6.54 


4.83 


0.97 


5.77 


1.60 


76.04 


0.65927 


0.37429 


NO 


0.008 


10.00 


5.70 


4.85 


0.97 


5.02 


1.60 


76.14 


0.75627 


0.42821 


NO 


0.008 


10.00 


4.99 


4.88 


0.97 


4.40 


1.60 


76.28 


0.86047 


0.48537 


NO 



be seen as the signature of a supernova explosion that changed all 
binary system parameters. Since PSR J 1 8 1 1 — 1736 is a recycled 
pulsar, the orbit must have been circularised during the recycling 
process (Battacharya & van den Heuvel 1991) and its orbital eccen- 
tricity could have been produced by a second supernova explosion, 
i.e. that of the companion star. In this case, both the spin period 
and the orbital eccentricity values of PSR J181 1 — 1736 would be 
the highest among DNSs and consistent with the correlation be- 
tween these two parameters observed in such systems (Faulkner et 
al. 2005) and recovered under the hypothesis of a low-amplitude 
neutron star kick (a v ~ 20km s~ ) at birth (Dewi et al. 2005). 

We investigated whether the MS companion scenario would 
be indeed compatible with the pulsar spin and orbital parameters. 
We compared in Fig. 5 the eccentricity, spin period, and orbital 
period of PSR J181 1 — 1736 to those of binary pulsars with iden- 
tified companions, whose masses are in the same range as the 
PSR J 1 8 1 1 — 1736 companion. We selected our sample from the 
Australia Telescope National Facility (ATNF) pulsar data bas^EI 
(Manchester et al. 2005). Our sample is summarised in Table 3. In 
our analysis, we focused on the comparison with recycled binary 
pulsars only, whose evolutionary path can be compared to that of 
PSR J181 1 — 1736. As seen, the only known pulsar-MS star sys- 
tem in the selected companion mass range is PSR J1903+0327 
(Khargharia et al. 2012). However, this is a fully-recycled ms- 
pulsar (P s = 2.5 ms), whereas PSR J1811-1736 is a mildly- 

5 http://www.atnf.csiro.au/people/pulsar/psrcat/ 



recycled pulsars with a much longer spin period (P s = 104 ms) 
and a much shorter orbital period (P or b = 18.8 d) than PSR 
J1903+0327 (P orb = 95.2 d). Moreover, the eccentricity of PSR 
J181 1-1736 (e = 0.828) is much larger than PSR J1903+0327 
(e = 0.436). Thus, there are no known pulsar-MS star systems in 
the selected companion mass range with spin and orbital param- 
eters comparable to those of PSR J 1 8 1 1 — 1736. This might sug- 
gest that such systems, if they do exists, are quite rare, although 
the very small sample currently available prevents us to draw firm 
conclusions. There is one pulsar system, PSR J0514— 4002A in 
the globular cluster NGC 1851 (Freire et al. 2004), with a possible 
white dwarf (WD) companion (Freire et al. 2007), that has both 
orbital period and eccentricity comparable to PSR J181 1 — 1736 
(Fig. 5, lower panel). However, like PSR J1903+0327, also PSR 
J0514— 4002A is a fully-recycled pulsar, with a spin period P s — 
4.99 ms. Moreover, since PSR J05 14— 4002A is in a globular clus- 
ter, its original companion might have been exchanged through a 
close encounter with another star in the cluster. Thus, the evolu- 
tionary history of this system might not be directly comparable to 
PSR J1811-1736. Afirm classification of the PSR J0514-4002A 
companion would help to evaluate the pulsar-MS star scenario for 
PSR J1811-1736. 



3.3 Constraints on the nature of the companion star 

If the companion star of PSR J181 1 — 1736 is undetected in the 
UKIDSS data, the derived upper limits on its flux can be used to 
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Table 2. Same as Tab. [T] but for ages smaller than 10 Gyrs. 
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0.85842 


0.60814 


NO 


0.020 


9.98 


5.70 


4.76 


1.06 


5.04 


1.51 


63.27 


0.98276 
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NO 
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NO 
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0.49420 


NO 


0.008 
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5.70 
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1.56 
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0.86716 


0.56498 


NO 


0.008 
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4.99 
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1.01 


4.41 


1.56 


69.62 


0.98515 


0.64185 


NO 


0.020 


9.92 


6.54 


4.74 


1.10 


5.77 


1.47 


59.57 


0.92229 


0.67548 


NO 


0.020 


9.92 


5.70 


4.77 


1.10 


5.03 


1.47 


59.57 


1.05797 


0.77486 


NO 


0.020 


9.92 


4.99 
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0.88167 


NO 


0.008 
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NO 


0.008 


9.92 
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1.02 


5.03 
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0.89589 


0.59975 


NO 
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4.99 
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0.68222 


NO 


0.020 
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1.46 
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0.69644 


NO 


0.020 
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4.79 
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NO 


0.020 
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4.99 
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58.53 
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NO 
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NO 
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1.03 


5.03 
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66.30 


0.92585 


0.63339 


NO 
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NO 


0.020 


9.80 


6.54 


4.49 


1.18 


5.24 


1.39 


53.51 


1.13065 


0.86080 


NO 
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NO 
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NO 


0.020 
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4.99 
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YES 


0.020 


9.60 


6.54 


4.84 
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1.23 


44.99 


1.15094 


0.90651 


NO 


0.020 


9.60 


5.70 


4.84 


1.34 


5.05 


1.23 


45.09 


1.33147 


1.04838 


YES 


0.020 


9.60 


4.99 


4.90 


1.33 


4.48 


1.24 


45.22 


1.49824 


1.17923 


YES 



constrain its nature. From the interstellar extinction coefficients of 
Fitzpatrick (1999), an E(B - V) = 5.7 ± 1.9 would correspond 
to Aj ~ 3.3-6.6, Ah ~ 2.0^1.0 and A K ~ 1.4-2.8. From our 
derived detection limits (J~ 20.5, H~ 19.4 and K~ 18.6) these 
values imply extinction-corrected fluxes of Jo > 13.9, Ho > 15.4, 
and A'o > 15.8, where we conservatively assumed the largest val- 
ues of the interstellar extinction. At the estimated pulsar distance 
(D = 5.7^0 71 kpc) these values correspond to absolute magni- 
tudes Mj > 6.9, M H > 8.4, and M K > 8.8, allowing us to rule 
out a companion of spectral type earlier than a mid-to-late M-type 
MS star. Thus, our constraints on the companion star are far more 
compelling that those derived by Mignani (2000) on the basis of 
the DSS data alone. Moreover, those constraints should be now re- 
vised upward since the reddening towards the pulsar measured in 
this work is at least twice as large than assumed by Mignani (2000) 
from the Galactic extinction maps (Hakkila et al. 1997). As a mat- 
ter of fact, these were the only resources available at the time to 
determine the reddening in the pulsar direction. 2MASS data of 
the pulsar field, which could be used to determine the reddening 



from the CMD technique, as we did with the UKIDSS data, were 
only released after the Mignani (2000) paper was published. Our 
new limits are also not deep enough to rule out a WD compan- 
ion star. As done in the previous section, we investigated whether 
the PSR J 18 1 1 — 1736 spin and orbital parameters would fit those 
of recycled pulsar-WD systems. As seen from Fig. 5, most recy- 
cled pulsar-WD systems have circular orbits and both orbital and 
spin periods shorter than PSR J181 1 — 1736. There is only one recy- 
cled pulsar-WD system, PSR J1750-3703A in the globular cluster 
NGC 6397 (D'Amico et al. 2001), that has both spin (P s = 0.111 
s) and orbital parameters (P rb = 17.33 d; e = 0.712) close to 
PSR J1811-1736. Thus, it is possible that PSR J1811-1736 is, 
indeed, a pulsar-WD system and not a DNS. However, since also 
PSR J1750— 3703A is in a globular cluster, the same caveats as dis- 
cussed above for PSR J0514— 4002A apply in this case. The iden- 
tification of other pulsar-WD systems with spin and orbital param- 
eter close to PSR J 1 8 1 1 — 1736, but located in the Galactic plane, 
would help to determine the nature of the PSR J1750— 3703A sys- 
tem. We note that the other pulsar-WD system PSR B 2303+46 
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Table 3. Name, timing (P B , P s , tsd, B Buli ) and orbital parameters (P or b, e) for binary pulsars with identified companions and masses in the same range 
as the PSR J1811 — 1736 companion. The sample has been selected from the ATNF pulsar data base (Manchester et al. 2005). Pulsar names are sorted in 
right ascension. The values of the companion masses Mq are either directly measured or calculated from other parameters (e.g. post-keplerian parameters). 
The two extremes of the mass range for Mc are computed assuming an inclination angle i = 90° and i = 60°, respectively, and a neutron star mass 
of 1.35 Mq. i = 60° is the median orbital inclination, for which there is an equal probability P of having an inclination smaller or larger than 60°, i.e. 
P(i < 60° ) = P(i > 60° ) = 0.5. The last two columns indicate the companion type (WD, NS, MS) and whether the pulsar system is ordinary or recycled. 
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Ps 


Ps 




T"SD 


^surf 


-forb 


e 




M c 


Companion 


Ordinary (O)/ 




(s) 


(ss- 1 


) 


(yrs) 


(G) 


(d) 






(Mq) 




Recycled (R) 


J0514— 4002A 


0.004991 


1.17 10" 


-21 


6.75 10 10 


7.73 10 7 


18.7852 


8.880 10" 


- 1 


0.90-1.11 


WD 


R 


B0655+64 


0.195671 


6.85 10" 


- 19 


4.52 10 9 


1.17 10 10 


1.0287 


7.500 10"" 


-6 


0.66-0.80 


WD 


R 


J0737-3039A 


0.022699 


1.76 10" 


- 18 


2.04 10 s 


6.40 10 9 


0.1023 


8.778 10" 


_2 


1.24890 


NS 


R 


J1022+1001 


0.016453 


4.33 10" 


-20 


6.01 10 9 


8.55 10 s 


7.8051 


9.700 10" 




1.05000 


WD 


R 


Jl 141— 6545 


0.393899 


4.31 10" 


-15 


1.45 10 6 


1.32 10 12 


0.1977 


1.719 10" 


-1 


1.02000 


WD 


O 


Jl 157—5112 


0.043589 


1.43 10" 


-19 


4.83 10 9 


2.53 10 9 


3.5074 


4.024 10" 


-4 


1.18-1.46 


WD 


R 


J1337-6423 


0.009423 


1.95 10" 


-19 


7.64 10 8 


1.37 10 9 


4.7853 


2.004 10" 


.5 


0.78-0.95 


WD 


R 


J1435-6100 


0.009348 


2.45 10" 


-20 


6.05 10 9 


4.84 10 s 


1.3549 


1.047 10" 


■5 


0.88-1.08 


WD 


R 


J1439-5501 


0.028635 


1.42 10" 


-19 


3.20 10 9 


2.04 10 9 


2.1179 


4.985 10" 


■5 


1.11-1.38 


WD 


R 


J1454-5846 


0.045249 


8.17 10" 


-19 


8.78 10 s 


6.15 10 9 


12.4231 


1.898 10" 


3 


0.86-1.05 


WD 


R 


J1518+4904 


0.040935 


2.72 10" 


-20 


2.39 10 1() 


1.07 10 9 


8.6340 


2.495 10" 


-1 


0.82-0.99 


NS 


R 


J1528-3146 


0.060822 


2.49 10" 


-19 


3.87 10 9 


3.94 10 9 


3.1803 


2.130 10" 


-4 


0.94-1.15 


WD 


R 


B1534+12 


0.037904 


2.42 10" 


-18 


2.48 10* 


9.70 10 9 


0.4207 


2.737 10" 


-1 


1.35000 


NS 


R 


J1750-3703A 


0.111601 


5.66 10" 


-18 


3.12 10 s 


2.54 10 1() 


17.3343 


7.124 10" 


-1 


0.58-0.69 


WD 


R 


J1756-2251 


0.028462 


1.02 10" 


-18 


4.43 10 s 


5.44 10 9 


0.3196 


1.806 10" 


-1 


1.10-1.35 


NS 


R 


J1802-2124 


0.012648 


7.26 10" 


-20 


2.76 10 9 


9.69 10 s 


0.6989 


2.474 10" 


■6 


0.78000 


WD 


R 


J1807-2459B 


0.004186 


8.23 10" 


-20 


8.06 10 s 


5.94 10 s 


9.9567 


7.470 10" 


-1 


1.20640 


WD 


R 


J1811 — 1736 


0.104182 


9.01 10" 


-19 


1.83 10 9 


9.80 10 9 


18.7792 


8.280 10" 


-1 


0.85-1.04 


NS 


R 


B1820-11 


0.279829 


1.38 10" 


-15 


3.22 10 6 


6.29 10 11 


357.7620 


7.946 10" 


-1 


0.65-0.78 


WD 


O 


J 1829+2456 


0.041010 


5.25 10" 


-20 


1.24 10 1() 


1.48 10 9 


1.1760 


1.391 10" 


-1 


1.26-1.57 


NS 


R 


J1903+0327 


0.002150 


1.88 10" 


-20 


1.81 10 9 


2.04 10 8 


95.1741 


4.367 10" 


-1 


1.03000 


MS 


R 


J 1906+0746 


0.144072 


2.03 10" 


-14 


1.13 10° 


1.73 10 12 


0.1660 


8.530 10" 


-2 


0.80-0.98 


NS 


O 


B1913+16 


0.059030 


8.63 10" 


-18 


1.08 10 8 


2.28 10 1() 


0.3230 


6.171 10" 


-1 


1.3886 a 


NS 


R 


B2127+11C 


0.030529 


4.99 10" 


-18 


9.70 10 7 


1.25 10 10 


0.3353 


6.814 10" 


-1 


1.354 6 


NS 


R 


B2303+46 


1.066371 


5.69 10" 


-16 


2.97 10 7 


7.88 10 11 


12.3395 


6.584 10" 


-1 


1.16-1.43 


WD 


O 



a ' b The values of the companion mass is taken from Weisberg et al. (2010) and Jacoby et al. (2006), respectively, and are not yet implemented in the ANTF 

pulsar data base 



(Dewey et al. 1985; Thorsett et al. 1993; van Kerkwijk & Kulka- 
rni 1999) with orbital parameters (P orb = 12.33 d; e = 0.658) 
similar to PSR J181 1 — 1736 (Fig. 5, lower panel) is not a recycled 
pulsar, which means that it is either on a different evolutionary path 
or evolutionary stage. 



4 SUMMARY AND CONCLUSIONS 

Using UKIDSS near-IR images, we detected a star, with magni- 
tudes J = 18.61±0.37, H = 16.65±0.03, and K = 15.46±0.02, 
within the 3a radio position uncertainty of PSR J181 1 — 1736. In 
order to determine the star's spectral type, we estimated the red- 
dening along the line of sight from the comparison of the CMDs 
of the stellar field to those of the Baade's Window, also built us- 
ing UKIDSS data. The reddening turns out to be at least twice 
as large as expected from the Galactic extinction maps. At a pul- 
sar distance of ~ 5.5 kpc, and for the estimated reddening of 
E(B — V) ~ 4.9, the star detected near to the radio position 
could be either a MS star close to the turn-off or a lower RGB 
star. The inferred mass (w 1.3M Q ) and radius (w 5R@) of this star 
could be compatible with the pulsar mass function, the constraints 
on the pulsar and companion masses, and the lack of radio eclipses 
near superior conjunction. Thus, it is possible that this star is the 
pulsar companion, which would reject the DNS scenario for PSR 



Jl 8 1 1 — 1736. if this is the case, this might be the first known ex- 
ample of a mildly-recycled pulsar-MS star system with companion 
mass in the ~ 0.9-1. 5M range, high eccentricity (e > 0.8), and 
spin and orbital periods in the explored range. However, we note 
that the computed chance coincidence probability of the candidate 
companion with the proper motion-corrected 3cr radio position is 
P ~ 0.27, which suggests that it might, instead, be an unrelated 
field star. A conclusive piece of evidence to prove/disprove the as- 
sociation would come from IR spectroscopic observations of the 
candidate companion star along the orbital phase of the binary sys- 
tem and the comparison of its velocity curve with that predicted 
by the orbital parameters of PSR J181 1 — 1736. Were the star con- 
firmed to be its companion, this would drive new theoretical studies 
on the birth and evolution of neutron stars in binary systems, the 
formation of mildly-recycled radio pulsars, and the amplitude of 
neutron star kicks imparted by supernova explosions. On the other 
hand, were the star proved to be an unrelated field star, the iden- 
tification of PSR J1811 — 1736 as a DNS would remain an open 
issue. Our near-IR detection limits with UKIDSS only rule out a 
companion of spectral type earlier than a mid-to-late M-type MS 
star. Deeper observations with 8m-class telescopes would enable 
us to push these limits down by about 4 magnitudes in each band. 
This would still not be enough to rule out any possible companion 
other than a neutron star, though, with a WD still compatible with 
the deepest achievable limits, unless the reddening is much lower 
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Figure 5. Orbital eccentricity e vs. spin P B (top panel) and orbital period P or b (lower panel) for the binary pulsars in Table 3. Different symbols and 
colours correspond to different pulsar systems: DNSs (green diamonds); pulsar-WD (blue triangles); pulsar-MS (magenta filled circle). Squares highlight 
non-recycled pulsar systems. The location of PSR J181 1 — 1736 is marked by the circled red diamond. The inset shows a zoom of the plots around the location 
of PSR J 1 8 1 1 — 1736. Pulsar-systems falling closest to PSR J 1 8 1 1 — 1 736 in the e-P s and e-P or b planes are labeled and colour-coded. 



than estimated in the current work. However, as shown in Sectn. 
3.3, most recycled pulsar-WD systems in the Galactic plane do not 
fit the spin and orbital parameters of PSR J181 1 — 1736. The de- 
tection of PSR J1811 1736 in X-rays, like the other DNS PSR 
J1537+1 155 (Durant et al. 201 1), would be useful to independently 
constrain the reddening along the line of sight and put tighter con- 
straints on the absolute luminosity of the companion. Finally, new 
radio observations of PSR J 1 8 1 1 — 1736 would be important to de- 
rive an updated radio-timing position and precisely measure the 
pulsar proper motion, for which we could only obtain a 2a mea- 
surement using the current radio observation data base. A more pre- 
cise radio- timing position of PSR J181 1 — 1736 will, then, enable 
us to revisit its association with its candidate companion star. 
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